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ABSTRACT 

The Advanced Camera for Surveys on board the Hubble Space Telescope has been used to obtain deep 
high-resolution images of the giant elliptical galaxy NGC 3610, a well-established remnant of a dissipative 
galaxy merger These observations supersede previous, shallower observations which revealed the presence 
of a population of metal-rich globular clusters of intermediate age ('^1.5-4 Gyr). We detect a total of 580 
cluster candidates, 46% more than from the previous WFPC2 images. The new photometry strengthens the 
significance of the previously found bimodality of the color distribution of clusters. Peak colors in V—I are 0.93 
± 0.01 and 1.09 ± 0.01 for the blue and red subpopulations, respectively. The luminosity function of the inner 
50% of the metal-rich ('red') population of clusters differs markedly from that of the outer 50%. In particular, 
the luminosity function of the inner 50% of the red clusters shows a flattening consistent with a turnover that 
is about 1.0 mag fainter than the turnover of the blue clusters. This is consistent with predictions of recent 
models of cluster disruption for the age range mentioned above and for metallicities that are consistent with the 
peak color of the red clusters as predicted by population synthesis models. The radial surface density profile 
of red clusters follows that of the underlying galaxy light more closely than in 'normal' elliptical galaxies, 
which is consistent with the intermediate-age nature of the red clusters. We determine the specific frequency 
of clusters in NGC 3610 and find a present-day value of Sn = 1-4 ± 0.6. Using published age estimates for 
the diffuse light of NGC 3610 as well as cluster disruption models, we estimate that this value will increase to 
Sn = 3.8 ± 1.7 at an age of 10 Gyr, which is consistent with typical Sn values for 'normal' elliptical galaxies. 
Our findings constitute further evidence in support of the notion that metal-rich cluster populations formed 
during major mergers involving gas-rich galaxies can evolve dynamically (through disruption processes) into 
the red, metal-rich cluster populations that are ubiquitous in 'normal' giant ellipticals. 

Subject headings: galaxies: star clusters — galaxies: elliptical and lenticular, cD — galaxies: individual 
(NGC 3610) — galaxies: interactions 



1. INTRODUCTION 

Recent deep imaging studies of 'normal' giant elliptical 
galaxies with the Hubble Space Telescope (HST) and large, 
ground-based telescopes have shown that these galaxies usu- 
ally contain rich globular cluster (GC) systems with bimodal 
color distributions (e.g., Zepf & Ashman 1993; Kundu & 
Whitmore 2001; Larsen et al. 2001; Peng et al. 2006). Typi- 
cally, roughly half of the GCs are blue, and half red. Follow- 
up spectroscopy with 8-m class telescopes has revealed that 
both 'blue' and 'red' GC subpopulations are typically old 
(> 8 Gyr, Forbes et al. 2001; Cohen, Blakeslee, & Cote 2003; 
Puzia et al. 2005) implying that the bimodahty is mainly due 
to differences in metallicity. The colors of the 'blue' GCs 
are usually similar to those of metal-poor halo GCs in the 
Milky Way and M31, while the mean colors of the 'red' GCs 
are similar to those of the diffuse light of their host galaxies 
(e.g., Geisler, Lee, & Kim 1996; Forbes, Brodie, & Grillmair 
1997). Hence, the nature of the 'red', metal-rich GCs is likely 
to hold important clues to the star formation history of their 

' Based on observations with the NASA/ESA Hubble Space Telescope, ob- 
tained at the Space Telescope Science Institute, which is operated by AURA, 
Inc., under NASA contract NAS5-26555. 



host galaxies. 

One environment known to produce metal-rich GCs and 
bimodal color distributions is that of vigorous starbursts in- 
duced by mergers of gas-rich galaxies. Massive young GCs 
have been commonly found in mergers and young merger 
remnants using HST observations (e.g., Holtzman et al. 1992; 
Schweizer 2002, and references therein). Follow-up spec- 
troscopy has confirmed the ages (and in three cases even the 
high masses, Maraston et al. 2004; Bastian et al. 2006) of 
these young clusters predicted from their colors and luminosi- 
ties (e.g., Schweizer & Seitzer 1998). Their metallicities tend 
to be near solar, as expected for clusters formed out of en- 
riched gas in spiral disks. A natural interpretation of these 
data is that the metal-rich GCs in 'normal', old giant ellipti- 
cals formed in gas-rich mergers at 2; > 1, and that the forma- 
tion process of giant ellipticals with significant populations of 
metal-rich GCs was similar to that in galaxy mergers observed 
today (e.g., Schweizer 1987; Ashman & Zepf 1992). 

An often quoted hurdle for this 'formation by merger' sce- 
nario has been the marked difference in the luminosity func- 
tions (hereafter LFs) of old vs. young GC systems (e.g., van 
den Bergh 1995). The LFs of old GC systems of 'normal' 
galaxies are well fit by Gaussians peaking at My ~ —7.2 
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mag with a dispersion of cr ~ 1.3 mag (e.g., Whitmore 1997), 
while young GC systems in mergers and young remnants have 
power-law LFs with indices of a ~ — 2 (e.g., Whitmore et al. 
1999). Intermediate-age merger remnants with ages of 1 - 5 
Gyr are ideal probes for studying the long-term dynamical ef- 
fects on GC systems formed during a major merger. Such 
galaxies are still identifiable as merger remnants through their 
morphological fine structure (e.g., Schweizer & Seitzer 1992), 
yet they are old enough to ensure that substantial dynamical 
evolution of the globular clusters has already occurred. Re- 
cent calculations of dynamical evolution of GCs (including 
two-body relaxation, tidal shocking, and stellar mass loss) 
show that the least massive GCs are expected to disrupt first as 
galaxies age, which can gradually transform the initial power- 
law LFs into LFs with Gaussian-like peaks or turnovers (e.g., 
Vesperini 1998; Baumgardt 1998; Fall & Zhang 2001; Ves- 
perini & Zepf 2003; Prieto & Gnedin 2006). 

Recent HST studies of candidate intermediate-age merger 
remnants have revealed that their 'red' GC subpopulations 
show LFs consistent with power laws (as expected if formed 
during a recent merger event; Goudfrooij et al. 2001b; Whit- 
more et al. 2002, hereafter W02). However, the combina- 
tion of field of view and sensitivity of the HST instrumen- 
tation was not quite efficient enough to allow a detection of a 
tumover in the LFs until just a few years ago. The first empir- 
ical evidence for dynamical evolution of red GC LFs from a 
power law toward a Gaussian-like shape in an intermediate- 
age early-type galaxy was facilitated by the unprecedented 
sensitivity of the Advanced Camera for Surveys (ACS) aboard 
HST. Based on ACS observations, Goudfrooij et al. (2004) 
found that the outer half of the red GC system of the 3-Gyr- 
old merger remnant NGC 1316 still shows a power-law LF, 
whereas the inner half shows a LF with a turnover character- 
istic of old GC systems. 

In an effort to test how common this result may be among 
nearby intermediate-age ellipticals, the current paper reports 
on new ACS imaging of NGC 3610. Several optical imag- 
ing studies have shown that NGC 3610 is a bona fide inter- 
mediate-age merger remnant. Its colors are anomalously blue 
for a giant elliptical galaxy of its luminosity (Schweizer & 
Seitzer 1992; Goudfrooij et al. 1994a; Idiart, Michard, & de 
Freitas Pacheco 2002), suggesting recent star formation ac- 
tivity. Its inner region features a small, strongly twisted inner 
disk (Scorza & Bender 1990; Whitmore et al. 1997). Such 
disks are a natural byproduct of gas-rich galaxy mergers in 
simulations (Barnes 2002). The galaxy's outer envelope is 
extremely boxy relative to its inner body and shows several 
non-concentric plumes and tails (Seitzer & Schweizer 1990) 
which are most likely remnants of tidal perturbations. In fact, 
NGC 3610 has the highest "fine structure index" (indicative 
of structures believed to be resulting from galactic mergers) 
of all 74 early-type galaxies in the sample of Schweizer & 
Seitzer (1992). Finally, its spheroid is well fit by a de Vau- 
couleurs' r^/^ profile and has a surprisingly small effective 
radius for the galaxy's luminosity (13" or 2.1 kpc in the B 
band, Burstein et al. 1987). Overall, these features are consis- 
tent with NGC 3610 being the product of a dissipative merger 
with incomplete dynamical relaxation. 

In this paper we adopt a distance of 33.9 Mpc for 
NGC 3610, corresponding to a distance modulus of (m — 
M)o = 32.65 mag as measured from surface brightness fluc- 
tuations (Tonry et al. 2001). The Galactic foreground red- 
dening towards NGC 3610 is >ly = 0.00 and 0.03 according 



to Burstein & Heiles (1984) and Schlegel et al. (1998), re- 
spectively. In accordance with most literature on photometry 
of NGC 3610 (e.g., Goudfrooij et al. 1994a; Whitmore et al. 
1997, 2002; Michard 2005), we adopt the Burstein & Heiles 
(1984) value in this paper. We do not make any corrections 
for reddening internal to NGC 3610, given the absence of 
any patchy absorption or color changes in optical images of 
NGC 3610 (Goudfrooij et al. 1994b; Whitmore et al. 1997) 
and the non-detection of far-infrared emission (Knapp et al. 
1989). 

2. OBSERVATIONS 

NGC 3610 was observed with HST on June 29th, 2003, 
using the wide-field channel (WFC) of ACS and the F435W, 
F555W, and F814W filters, with total exposure times of 1080 
s, 6410 s, and 6060 s, respectively. The F555W and F814W 
exposure times were chosen to reach a similar depth for an 
object of V—I = 1.0, the mean V—I color of globular cluster 
systems of giant elliptical galaxies (e.g., Kundu & Whitmore 
2001; Larsen et al. 2001). The F435W exposures were mainly 
included to provide a large color baseline to allow accurate 
age/metallicity measurements for the bright GCs, which does 
not require reaching very faint levels. Hence, the F435W 
exposures were relatively short and taken at only two posi- 
tions ("dithers") offset from each other by 3.5 columns of the 
WFC CCD. In contrast, the final F555W and F814W images 
were constructed from long (19-26 min) exposures taken at 
5 dither positions, supplemented by one short exposure taken 
in each filter to avoid saturation of the bright central regions. 
The individual images in each band were combined using the 
task MULTIDRIZZLE within PyRAF/stsdas v3.3, with task 
parameters set to pixfrac = 0.8 and scale = 1.0. Saturated pix- 
els in the long exposures were replaced by those in the short 
exposures while running MULTIDRIZZLE by setting the ap- 
propriate data quality flag in the affected pixels of the long 
exposures. The gain in sensitivity, spatial resolution, and spa- 
tial coverage of the new ACS images relative to the available 
WFPC2 images aUowed us to detect roughly 50% more GC 
candidates than the WFPC2 images used in W02 did. 

3. CLUSTER SELECTION AND PHOTOMETRY 
3.1. Cluster Selection 

Prior to performing source selection and photometry, the 
strongly varying light distribution of the galaxy was fitted and 
removed in order to minimize errors in the photometry due 
to any particular choice of object aperture and sky annulus 
(see below). Experimentation with various fitting methods 
revealed that any model fit to the galaxy surface brightness 
distribution is strongly hampered by the presence of a very 
prominent stellar disk in the central regions (e.g., Scorza & 
Bender 1990; Goudfrooij et al. 1994a; Whitmore et al. 1997, 
2002). Best results were achieved by approximating the dif- 
fuse light of the galaxy by applying a median filter with a 
41x41 pixeP kernel to the image. 

The first step in the selection of candidate GC targets was 
the application of the daofind task from the daophot-ii 
package (Stetson 1987) to an image prepared by dividing the 
sum of the F555W and F814W images (i.e., the V+I image) 
by the square root of the corresponding median-filtered im- 
age (thus yielding uniform shot noise characteristics over the 
whole image). The detection threshold was set at 3 cr above 
the background. The region around the central stellar disk was 
masked out in this step. To allow the detection of candidate 
GCs in the region occupied by the stellar disk, daofind was 
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applied again on the V + 1 image divided by the square root 
of a median-filtered version of it, but this time using a much 
smaller kernel (7x7 pixeP). From this last run of DAOFIND, 
we only included objects detected in the central area that was 
masked out in the first run into the final object list. 

At this stage, a total of 1472 objects were detected. Task 
PHOT was then applied to obtain fixed-aperture photometry 
for all those objects in the median-subtracted B, V, I, and 

V + 1 images, employing several aperture radii. Due to the 
substantial structure in the median-filtered images in the in- 
nermost region, aperture photometry for the objects in the in- 
nermost region was obtained from the direct images, employ- 
ing small apertures (radii for object/inner background/outer 
background were set at 2/3/5 pixels in this case; the photome- 
try of these inner objects is shown with open symbols in Fig- 
ure 4). 

Since the object list at this point still includes Galactic 
foreground stars, background galaxies, and any residual cos- 
mic rays and/or hot pixels in addition to the star clusters we 
are interested in studying, we devised an automated selec- 
tion procedure based on two size-related parameters and the 

V — I color. Quantitative size selection criteria were devel- 
oped based on globular clusters known to be associated with 
NGC 3610 by means of radial velocities (Strader, Brodie, & 
Forbes 2004), obvious bright foreground stars with visible 
diffraction spikes, and visual inspection of objects (especially 
near the extrema of the selection criteria). The first size pa- 
rameter used is a 'concentration index', defined as the differ- 
ence in F555 W magnitudes measured using aperture radii of 2 
and 5 pixels, hereafter denoted as F555W2-5. Visual inspec- 
tion of the objects remaining after restricting F555W2-5 to 
the range 0.29 - 0.83 mag showed that many extended back- 
ground galaxies, all obvious hot pixels, and some bright fore- 
ground stars were filtered out. However, this selection still 
left a substantial number of background galaxies in the tar- 
get list. The second size parameter used for the classification 
is the FWHM of the radial profile of the object. Rather than 
assuming a specific model of the radial profile (e.g., Gaus- 
sian or Moffat), we used a more robust, direct measure of the 
FWHM by computing the radial derivative of the azimuthally 
averaged enclosed flux profile F{R): 

F{R) = / P{r)rdr 
Jo 

where P{r) is the radial profile of intensity per unit area. The 
peak of dF{R)/ Rdr is found and the FWHM is twice the 
radius of the profile at half the peak value. As illustrated in 
Figure 1, we selected objects with 0.29 < F555W2-5 < 0.83 
mag and 1.75 < FWHM (pixels) < 3.50 in this step. 

At the distance of NGC 3610, typical GCs are only 
marginally resolved at the resolution of the drizzled ACS im- 
ages. This makes it somewhat difficult to distinguish fore- 
ground stars from GC candidates based on size parameters 
alone, especially at faint magnitudes where the absence of 
diffraction spikes hampers the visual identification of objects 
as foreground stars. Hence we use the V — I color as another 
discriminant between stars and GC candidates by applying a 
color selection cut of 0.5 < {V — I)o < 1.5. This interval 
constitutes the full range in V~I expected for GCs older than 
1 Gyr and —2.5 < [Z/H] < 0.5 according to the simple stel- 
lar population (SSP) models of Bruzual & Chariot (2003) and 
Maraston (2005), plus 0.2 mag of padding on the blue and red 
ends to allow for reasonable photometric errors. These con- 
centration, size, and color cuts yielded a list of 580 GC can- 
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Fig. 1. — Illustration of the size selection procedure. The FWHM of the 
object profile measured as discussed in Sect. 3.1 is plotted versus F555W2-5, 
the difference in F555W magnitudes measured using aperture radii of 2 and 
5 pixels. The region within dashed lines represents the parameter space used 
to select GC candidates in NGC 3610. 



didates. Any further contamination of this list by foreground 
stars and compact background galaxies will be addressed by 
statistical subtraction (see below). 

3.2. Photometry 

Aperture corrections were —0.26 (±0.01) mag in F435W 
(B), -0.24 (±0.01) in F555W {V) and -0.29 (±0.01) mag 
in F814W (/) for an aperture radius of 3 WFC pixels (i.e., 
0'.'15). These values were determined through measure- 
ments of several isolated, well-detected GC candidates lo- 
cated throughout the field. These GC candidates spanned the 
range of size parameters mentioned in the previous section. 
We considered correcting for the effect of imperfect charge 
transfer efficiency (CTE). However, the effect was found to be 
negligible (< 0.01 mag) according to the formulae of Riess & 
Mack (2004), even at the corners of the F435W image which 
have the lowest background level by far (^ 30 e^/pixel). A 
correction for CTE loss was therefore not applied. In order to 
facilitate comparisons of our photometry with predictions of 
various SSP models, we transformed the instrumental magni- 
tudes to Johnson/Cousins B, V and / using the calibrations 
derived by the ACS Instrument Definition Team (Sirianni et 
al. 2005). Photometry and astrometry for the 50 brightest GC 
candidates in the ACS frames are listed in Table 1 . 

Considering the sky coverage common to this study and that 
of W02, a comparison of the photometry in Table 1 with that 
of the brightest 50 GC candidates in the WFPC2 data of W02 
(i.e., their Table 1) shows a weighted mean difference (in the 
sense "our minus W02") in V magnitudes of —0.01 with an 
rms scatter of 0.12 mag. (The weighting was done by means 
of the inverse variance.) The difference in weighted mean 
V — I colors is —0.02 with a rms scatter of 0.11 mag. This 
comparison excludes objects #1,4, and 8 of W02, which are 
all located within the central disk of NGC 3610. A careful 
check of the ACS images does not show any evidence for the 
presence of objects on top of the bright disk light at these 
locations. 

Completeness corrections were performed using the V + I 
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TABLE 1 

Photometry AND astrometryof the 50 brightest GC candidates in NGC 3610 from the ACS data. The object list is sorted on V 

MAGNITUDE (BRIGHTEST FIRST). 
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image, i.e., the one used to determine the initial Ust of GC can- 
didates, using routines in the DAOPHOT package. We added 
artificial objects (in batches of 100) with a radial intensity pro- 
file derived from a fit to real GC candidates in the frame to the 
V + 1 image, for five different background levels and several 
magnitude intervals. The color of the objects was set equal to 
the median color of GC candidates. Prior to running DAOFIND 
to determine the fraction of retrieved artificial objects, the im- 
age was divided by the square root of its median-filtered ver- 



sion, as was done during the initial object finding procedure. 
Artificial objects retrieved by daofind were also subjected to 
the size selection criteria mentioned above using the median- 
filter-subtracted V+I image before including them in the ob- 
ject counts^. Figure 2 shows the resulting completeness func- 

^ To reflect the different method of object photometry in the innermost 
regions (see Section 3.1), the photometry of the artificial objects at the highest 
backgroimd level was done on the direct image instead of the median-filter- 
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tions. 

The two main reasons for using the V + I image to de- 
termine the initial source list and completeness corrections 
(rather than the V or / image, which is often done) are (i) 
the greater depth reached by the V+I image, and (ii) the fact 
that the transformation from instrumental magnitudes to John- 
son/Cousins V + 1 is significantly less dependent on source 
color than the transformations to either V or /, which helps 
during the completeness correction. The latter point is illus- 
trated in Figure 3 which shows the dependence of the pho- 
tometric zeropoints for V and V + 1 on the V — I color 
The zeropoints were derived using the SYNPHOT package 
within IRAF/stsdas. Synthetic spectra of stellar types rang- 
ing from F2V to MOIII from the Bruzual & Chai'lot (1993) 
library were convolved with Johnson/Cousins V and / fil- 
ter passbands to yield Vega-based magnitudes, and with the 
ACS/WFC F555W and F814W passbands to yield STMAG 
magnitudes. As Figure 3 shows, the dependence of the V + I 
zeropoint onV — I color is much smaller than that of the V 
zeropoint. In the context of the two common subpopulations 
of clusters in giant ellipticals, which typically have V — I ^ 
0.95 and 1 .2 for the blue and red subpopulations, respectively, 
the difference in V+I zeropoint is only 0.006 mag (vs. 0.016 
mag for the V zeropoint). Hence, the determination of com- 
pleteness fractions using V + I needs only be conducted for 
one object color^. 

4. PROPERTIES OF THE GLOBULAR CLUSTER SYSTEM 
4. L Color Distributions 

The top panels of Figure 4 show the B vs. B — I and V 
vs. V—I color-magnitude diagrams (hereafter CMDs) for GC 
candidates in NGC 3610. For reference, the region occupied 
by old halo GCs in our Galaxy (taken from the database of 
Harris 1996), when placed at the distance of NGC 3610, is 
indicated by dashed lines in both CMDs. 

The CMDs show a color distribution of GC candidates that 
is clearly bimodal down to i? ^ 26.5 and V ^ 26.5 where the 
typical color uncertainties are ±0.22 mag in _B — / and ±0.10 
mag in V — I. The 'blue' sub-population of GCs has col- 
subtracted image. 

^ Note however that these advantages for using the summed (V+I) image 
are only applicable if the individual V and / images are of similar depth for 
the median color of the objects of interest. 
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Fig. 3. — The relation of the photometric zeropoints of V (open squares) 
and V + I (filled squares) with the V — I color Note that the dependence 
of the V + I zeropoint on V — / color is significantly smaller than that of 
the V zeropoint. Dashed lines show the zeropoints for typical 'blue' globular 
clusters found in NGC 3610 (V — I ^ 0.9), and dash-dot lines do so for 
typical 'red' globular clusters in NGC 3610 (V — I ^ 1.1). See discussion in 
Sect. 3.2. 



ors and luminosities that are consistent with those of old halo 
GCs in our Galaxy, while the (larger) 'red' population has a 
mean color that is ^0.4 mag redder than that of the 'blue' 
sub-population. The mean colors of each subpopulation seem 
constant throughout the sampled range of luminosities, tak- 
ing the color uncertainties as a function of brightness into ac- 
count. There is no obvious evidence for the 'red slant' effect 
seen among 'blue' GCs in brightest cluster ellipticals, i.e., the 
effect that the mean color of 'blue' GCs gets progressively 
redder (indicating a higher metallicity) with increasing GC 
luminosity (Harris et al. 2006; Strader et al. 2006). The V vs. 
V — I CMD looks quite similar to the corresponding diagram 
in W02 (i.e., their Figure 3), except that the color bimodality 
seems more obvious in the current CMD. 

The B—I and V—I color distributions of the GC candidates 
are shown in the bottom panels of Figure 4, for two brightness 
cuts. Prior to plotting, the color histograms of candidate GCs 
were corrected for contamination by compact background ob- 
jects. This correction was derived from the compact objects 
detected beyond a projected radius of 100" from the galaxy 
center The color histogram of the outer compact objects was 
smoothed with a five-bin kernel to diminish small-number 
noise and scaled to the area of the full ACS image prior to sub- 
traction from the color histogram of the full ACS image. The 
effect of this correction is shown in Figure 4. While the sur- 
face number density of GC candidates beyond 100" from the 
galaxy center seems to be close to an asymptotic background 
value (see Sect. 4.3 below), there may still be a (small) num- 
ber of true halo GCs included in the number density estimate 
of compact background galaxies. Hence the correction for 
compact background galaxies may be slightly overestimated, 
and the true color distribution of GCs in NGC 3610 may lie 
slightly above (but most likely very close to) the solid lines in 
the bottom panels of Figure 4. 

The unbiased probability of a bimodal color distribution (as 
opposed to a unimodal one) was estimated using the Kaye's 
Mixture Model (KMM) test (McLachlan & Basford 1988; 
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Fig. 4. — Top panels: B vs. B—I (left) and V vs. V—I (right) color-magnitude diagrams for GC candidates in NGC 3610. The GC candidates with projected 
galactocentric radii < 12" are shown as open circles, while the others are shown as filled circles. Regions within dashed lines represent the magnitude and color 
ranges for old halo GCs in the Milky Way, placed at the distance of NGC 3610. Representative error bars are shown on the right-hand side of the diagrams. 
Bottom panels: Distributions oi B — I (left) and V — I (right) colors for GC candidates in NGC 3610. The bottom left panel shows all GCs with B < 26.5 
(open histogram), and GCs with B < 25.5 (filled histogram). The dashed fines represent GCs with B < 26.5 before correction for contamination by compact 
background galaxies. The coiTesponding histograms in the bottom right panel are for V < 26.5 (open histogram and dashed lines) and for V < 25.0 (filled 
histogram). 



Ashman, Bird, & Zepf 1994), by comparing fits of one Gaus- 
sian vs. two Gaussians to the color distribution using the 
"expectation-maximization" (EM) method. To gain leverage 
in the presence of noisy data, we consider the homoscedastic 
case of the KMM test where the standard deviation a is the 
same for both Gaussians. Using the KMM test on various ar- 
tificial color distributions, we determine that a two-Gaussian 
fit is a better fit to the data than a one-Gaussian fit if the p- 
value of the two-Gaussian fit is less than 0.05. We ran the 
KMM test on several subsets of the color distributions of GC 
candidates, using two magnitude limits and two upper limits 
of the projected galactocentric radius. Table 2 lists the results 
of the KMM tests. In all cases the distributions are judged 
to be bimodal at the 96-100% probability level. The mean 
colors of the blue peak are S — / = 1.61 ± 0.01 and V — I = 
0.93 ± 0.01, and those of the red peak are B~I = 1.92 ± 0.01 
and y-/= 1.09 ±0.01. 

The dependence of the color distribution of GC candidates 
on galactocentric radius is illustrated in Figure 5, which shows 
V—I color vs. projected galactocentric radius for all GC can- 
didates with V < 26.5. In addition. Figure 5 shows average 
colors of GC candidates in bins of log (radius) for three sam- 
ples: (i) all GC candidates, (ii) GC candidates in the 'blue' 
subpopulation(0.75 < V-I < 1.00), and f/Z/j GC candidates 
in the 'red' subpopulation (1.00 < V — I < 1.30). The aver- 
age color of the full sample of GC candidates gets bluer with 
increasing projected radius, as found before for many 'nor- 
mal' ellipticals (e.g.. Ashman & Zepf 1998, and references 
therein). Fig. 5 also shows that the average colors of the blue 
and red subpopulations stay constant with increasing radius 
(see Fig. 5), and that these results are not affected by contam- 



ination by compact background galaxies. Hence, the overall 
radial color gradient is due to a more centrally concentrated 
radial distribution of the red GCs with respect to their blue 
counterparts. This effect (which is quantified for NGC 3610 
in Section 4.3 below) was predicted by the 'formation by dis- 
sipative merger' scenario of Ashman & Zepf (1992), and con- 
firmed by recent studies of giant ellipticals hosting populous 
GC systems with bimodal color distributions (e.g., Geisler et 
al. 1996; Goudfrooij et al. 2001b; Sikkema et al. 2006). 

The peak colors of the blue and red subpopulations men- 
tioned above are compared with popular SSP models in Fig- 
ure 6. We considered the models of Anders & Fritze-v. 
Alvensleben (2003, using the Salpeter initial mass function 
(IMF)), Bruzual & Chariot (2003, using the Chabrier IMF), 
and Maraston (2005, using the Kroupa IMF) for this purpose. 
As to the observed peak B—I and V—I colors of the blue GCs, 
Figure 6 shows that all three models fit these quite well at an 
age of 14 Gyr and [ZfH] = —1.7, which is consistent with 
the mean [Z/H] of old Galactic halo GCs (e.g., Harris 1996). 
Given that the luminosity function of the blue GC subsystem 
is also consistent with that of the Galactic halo GCs (see be- 
low), we consider the blue GC subsystem to consist of old, 
metal-poor GCs for all intents and purposes. As to the red GC 
subsystem. Figure 6 shows that its peak B — I and V — I col- 
ors are consistent with various age/metallicity combinations. 
These are discussed in turn below. 

One option to explain the colors of the red GCs is that they 
are 'old' (age > 10 Gyr) GCs with -1.0 < [Z/H] < -0.5 
(the exact best-fit value for [Z/H] depending on the preferred 
SSP model), i.e., GCs with properties similar to those of bulge 
GCs in our Galaxy (e.g., Puzia et al. 2002; Maraston et al. 
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TABLE 2 

Results of the Kaye's Mixture Model (KMM) tests to the color distributions 
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^Outer galactocentric radius of GC sample under consideration. 

''Number of GCs in sample. 

■^Standai'd deviation of Gaussian fit by KMM test. 

''Fraction of red GCs in color distiibution indicated by KMM test. 
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Fig. 5. — Top panel: V — I color versus logarithm of projected galactocen- 
tric radius for GC candidates in NGC 3610 with V < 26.5. The plus signs 
represent the individual GC candidates. The large open black circles rep- 
resent the mean colors of GC candidates using bins of 0.25 in Log (radius) 
(except for the innermost bin, which represents Log (radius) < 1). The large 
filled red and blue circles represent the mean colors of the red and blue GC 
subpopulations, respectively, using the same radial binning as for the black 
circles. En'or bars drawn for the large circles represent mean errors of the 
mean. Values mentioned beside each red and blue circle indicate the local 
surface number density of red and blue GC candidates (in arcsec"^), respec- 
tively. Bottom panel: V — I color distribution of 'contaminating' compact 
background galaxies as estimated from the outer regions of the image, ex- 
pressed in surface number density units (arcsec~^). See text in Section 4.1 
for further details. 



2003), which may have been introduced into NGC 3610 if 
(at least) one of the merger progenitor spiral galaxies hosted 
a significant bulge component. The existence of old GCs 
among the 'red' GC subsystem in NGC 3610 was shown by 



Strader et al. (2004), who determined ages and metallicities 
of 1 1 GCs in NGC 3610 from spectra obtained with the Keck 
I telescope, using spectroscopic Lick indices. Based on SSP 
model fits to the Lick indices, Strader et al. argued that 7 out 
of the 9 'red' GCs in their sample were best fit by old ages 
(> 10 Gyr) and metallicities similar to those of GCs in the 
Galactic bulge. This option will be considered further in in- 
terpreting the LF of red GCs in Section 4.2 below. 

Another option that is consistent with the colors of the red 
GCs is that they are of intermediate age (~ 1-4 Gyr) and 
metal rich ([Z/H] > 0.0). Support for this option is provided 
by several empirical results. First of all, the LF of the red GC 
population in NGC 3610 is known to be well-fit by a power 
law down to well beyond the turnover magnitude of known 
'old' GC systems in giant ellipticals (W02; see also Sect. 4.2.2 
below). Second, the peak color of the red GCs in NGC 3610 
is significantly bluer than that of red GCs in 'normal' ellipti- 
cal galaxies at the luminosity of NGC 3610. This is illustrated 
in Figure 7 which shows the relation between the peak V — I 
colors of the red GC subpopulation in nearby elliptical galax- 
ies with clear bimodal color distributions versus the galax- 
ies' absolute _B-band magnitude from the study of Larsen et 
al. (2001). Figure 7 clearly shows that these two properties 
correlate strongly with each other, which has been generally 
interpreted as a galaxy luminosity-metallicity relation simi- 
lar to that among elliptical galaxies themselves, suggesting a 
close connection between the physical processes responsible 
for the formation of the red GCs and that of their parent galax- 
ies (e.g., Brodie & Strader 2006, and references therein). As- 
suming the same physical relation holds for NGC 3610 (and 
for the 3-Gyr-old merger remnant NGC 1316: Goudfrooij et 
al. 2001b, 2004), the fact that the peak color of their red GCs 
are clearly offset to the blue from this relation (see Figure 7) 
must be due to a younger age rather than a lower metallicity. 
Finally, the Lick indices of the 2 innermost red GCs in the 
spectroscopic sample of Strader et al. (2004) were best fit by 
an age of ^ 1.6 Gyr and metallicity [Z/H] ~ -1-0.3, consistent 
with Lick-system analysis of spectra of the diffuse light of the 
galaxy out to beyond its effective radius (Howell et al. 2004). 
Figure 6 shows that this age/metallicity combination fits the 
observed B — I and V — I colors quite well. 
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Fig. 6. — Time evolution ot B — I (top panels) and V ~I (bottom panels) color indices according to single-burst stellar population models for ages older than 
1 Gyr. Panels (a) and (b) show Bruzual & Chariot (2003) models, panels (c) and (d) show Anders & Fritze-v. Alvensleben (2003) models, and panels (e) and (f) 
show Maraston (2005) models. Model curves are plotted in black for the following metallicities: 0.02 solar (short-dashed lines), 0.2 solar (long-dashed lines), 
1.0 solar (solid lines), and 2.5 solar (dot-long-dashed lines). The horizontal red and blue dotted lines indicate the peak colors as determined by the KMM test for 
the red and blue GC subpopulations, respectively (see Sect. 4.1). The horizontal long-short-dashed lines above and below the red and blue dotted lines indicate 
the ±1(T widths of the two Gaussians fitted to the bimodal color distribution. The vertical blue dot-short-dashed lines in each panel indicate the 14 Gyr model 
age. The vertical red dot-short-dashed lines indicate the best-fit ages to the red peak color for solar metalHcity (right fine) and for 2.5 solar metallicity (left line), 
separately for each SSP model. The exception is the Maraston (2005) model where a range of ages (2-3 Gyr) fits the colors well for solar metallicity, and 2.5 
solar metallicity does not fit the data well for ages > 1 Gyr. Hence only one red vertical line is drawn for that case. 
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Fig. 7. — Peak V—I color of the red GC subpopulation versus absolute B- 
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bar is shown at the top left of the figure. The dashed fine represents a linear 
least-squares fit to the open squares. The data for NGC 3610 and NGC 1316 
(data for the latter was taken from Goudfrooij et al. 2004) are shown by filled 
squares. See discussion in Sect. 4.1. 



4.2. Cluster Luminosity Functions 

The LFs of the blue and red GC subpopulations are shown 
in Figure 8, and were determined as follows. We first sepa- 
rated the GC candidates into 'blue' (0.75 < V-I < 0.98) 
and 'red' (1.03 < V-I < 1.30) subsamples. Note that in 
this case, we avoided the overlapping region (0.98 < V—I < 
1.03) to minimize the mixing of the two subsamples. A com- 
pleteness correction was then applied to each GC by dividing 
its count (of 1) by the completeness fraction corresponding to 
its brightness and background level. The latter fraction was 
calculated from the functions shown in Figure 2, using bilin- 
ear interpolation in log(background) space. The LFs were fur- 
ther corrected for contamination by foreground stars and com- 
pact background objects using the smoothed and scaled LF 
of the objects detected at galactocentric distances R > 100" 
with the same color and size criteria as the GC candidates at 
B < 100". The LFs of those outer objects were smoothed 
by a five-bin median filter to avoid noise due to small-number 
statistics and scaled to the area of the R < 100" region prior 
to subtraction. The observed LFs corrected for completeness 
are shown as dotted lines in Figure 8, whereas the solid lines 
show the LFs after correction for contamination by compact 
background objects. 

As already reported by W02, the LFs of the (full) blue and 
red GC subsystems in NGC 3610 differ significantly from 
each other. We discuss the LFs of the two subpopulations 
separately below. 
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Fig. 8. — V-band luminosity functions (LFs) of GC candidates in NGC 3610. Panel (a): LF of the full blue subsystem. Panel (b): LF of the full red 
subsystem. Panel (c): LF of the outer 50% of the red subsystem. Panel (d): LF of the inner 50% of the red subsystem. Dotted histograms mark LFs corrected 
for incompleteness, and solid histograms mark LFs corrected for incompleteness and compact background galaxies. Histograms are filled for magnitude bins 
brighter than the overall 50% completeness limit of the sample plotted, and open beyond it. The smooth dashed curve in panel (b) is a power-law fit to the LF. 
Vertical arrows in all panels indicate the typical turnover magnitude for 'old' GC systems (i.e., My = —7.2). Note the appearance of a flattening consistent 
with a turnover in the LF of the inner 50% of the red subsystem, which is not present in the outer 50%. The smooth curves in panels (a) and (d) represent 
predicted LFs based on a combination of the Fall & Zhang (2001) GC disruption models and the SSP models of Maraston (2005). Curves are drawn for the 
following populations: ages of 1.5 (long-dashed line), 3.0 (short-dashed line), 6.0 Gyr (short-dash-dotted line), and 12 Gyr (solid line in panel (a)). [ZfH] = — 1.7 
is assumed for the 12 Gyr curve, whereas [Z/H] = 0.0 is assumed for the three younger ages. The normahzation of all model curves in panels (a) and (d) was 
done by a least-squares fit to the underlying histograms. 



4.2.1. The Blue Subpopulation 

The shape of the LF of the blue GC subsystem is simi- 
lar to that typically reported for GCs in 'normal' giant el- 
lipticals, i.e., approximately Gaussian with a turnover mag- 
nitude My ^ —7.2 (e.g., the compilation of Whitmore 
[1997] found My = -7.21 ± 0.26). A Gaussian fit to 
the LF of the blue GC subsystem of NGC 3610 results in 
a peak of A/y = —7.21 ± 0.15 and a standard deviation 
a = 1.22 ± 0.18. While the latter is consistent with that 
found among 'normal' giant early-type galaxies (e.g., Kundu 
& Whitmore [2001] found values between 1.1-1.5 mag), it 
is significantly larger than that found for NGC 3610 by W02, 
who found a = 0.66 ± 0.19. We ascribe this difference to a 
larger field of view and a larger sample of GC candidates in 
the current study relative to W02. 

4.2.2. The Red Subpopulation 

The LF for the full red GC system is consistent with a 
power law down to the 50% completeness limit of the data. 
A weighted least-squares fit of (t){L)dL oc L^dL to the LF 
yields an exponent a — —1.78 ± 0.06, consistent with the 
value found by W02 (a ^ -1.78±0.05). While this confirms 
the nature of the bulk of the red GC population of NGC 3610 
as having formed during a recent star-forming event, it does 
not provide significant evidence to suggest an evolution of the 
LF from a power law to a Gaussian. 

However, disruption timescales of the main dynamical ef- 
fects thought to be responsible for the Gaussian shape of the 
GC LFs of 'normal', old galaxies through preferential deple- 
tion of low-mass GCs (Vesperini 2001; Fall & Zhang 2001; 



Baumgardt & Makino 2003) depend on galactocentric dis- 
tance. The stronger tidal field in the central regions yields 
stronger tidal shocks than in the outskirts, while also impos- 
ing a more stringent zero-energy surface on the GCs so that 
stars can be more easily removed from the GCs through two- 
body relaxation (see also Gnedin & Ostriker 1997; Gnedin 
1997). Hence, the appearance of a turnover in the GC LF can 
be expected to first occur in the inner parts of the GC system. 
This was first shown to be the case by Goudfrooij et al. (2004) 
for the 3-Gyr-old merger remnant NGC 1316, where the in- 
ner half of the red GC system showed a significant flattening 
of the GC LF beyond the turnover expected (and found) for 
the blue, metal-poor GC system. In contrast, the outer half of 
the red GC system still rose strongly all the way to the detec- 
tion limit. With the 580 GC candidates in our ACS images of 
NGC 3610, we are now in a position to test this finding for this 
intermediate-age merger remnant with high enough statistical 
significance as well. To do so, we divide the red GC system 
into two equal-size parts sorted by projected galactocentric ra- 
dius. For our sample, this division ends up at 46'.'7, or 4.8 kpc 
from the galaxy center The resulting LFs are shown in panels 
(c) and (d) of Figure 8. The LF of the outer half of the red GC 
system still rises all the way to the detection limit, whereas 
the LF of the inner half of the red GC system shows a clear 
flattening consistent with a turnover at ~ 1.0 mag fainter than 
that of the blue GC system. These properties are quite similar 
to those found for the red GC system in NGC 1316 (Goud- 
frooij et al. 2004). 

We compare these results with predictions of GC disruption 
models by Fall & Zhang (2001, hereafter FZOl) in Figure 8. 
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Fig. 9. — V-band LFs of the red GC candidates in NGC 3610. Both panels show the same LF as in Figure 8d. Panel (a): FZOl model curves are drawn 
for the following populations of solar metallicity (using the Maraston (20(35) SSP models): ages of 1.5 (long-dashed line), 3.0 (short-dashed line), and 6.0 Gyr 
(short-dash-dotted line). Panel (b): FZOl model curves are drawn for the following populations: Age = 1.5 Gyr and [Z/H] = +0.3 (long-dash-dotted line) and 
Age = 12 Gyr and [Z/H] = —0.7 (dotted Hne). The normalization of all aforementioned model curves was done by a least-squares fit to the underlying histogram. 
To illustrate the estimated contribution to the red GC LF by a population of GCs associated with the bulge of a giant spiral galaxy similar to the Milky Way, the 
solid line represents a FZOl model curve for Age = 12 Gyr and [Z/H] = —0.7 for a population of 50 red GCs (see discussion in Sect. 4.2.2). For reference, a 
power-law fit to the full red GC population (cf. Figure 8b) is shown as a long-dashed curve in panel (b). 



The FZOl models, whose initial GC velocity distribution in- 
volves a radial anisotropy distribution that is consistent with 
kinematics of halo stars in the solar neighborhood, showed 
that an initial power law or Schechter mass function evolves 
in 12 Gyr to a peaked mass distribution similar to that ob- 
served for the Milky Way GC system. To predict the shapes 
of luminosity functions as a function of age, we use the FZOl 
model that uses a power-law initial mass function in combi- 
nation with the SSP models of Maraston (2005). We assume 
[Z/H] = — 1.7 for the 12 Gyr old population (similar to the 
median metallicity of halo GCs in our Galaxy), while we use 
[Z/H] = 0.0 for populations of ages 1.5,3, and 6 Gyr Prior to 
plotting, the amplitude of all resulting model LFs were nor- 
malized to the observed data down to the 50% completeness 
limit by weighted least-square minimization. The individual 
RMS values of these fits to the histograms are listed in Ta- 
ble 3. The peak magnitudes of the 1.5, 3, and 6 Gyr model 
LFs of solar metallicity span a range of only 0.3 mag in V, 
due to the counteracting effects of disruption of (preferentially 
low-mass) GCs and luminosity fading of the surviving GCs 
(see also W02). Hence the LF shapes of these three mod- 
els are quite similar to one another. However, as discussed 
in Sect. 4.1, recent spectroscopic investigations indicated that 
the intermediate-age GC population is of age 1 .5 Gyr and 
has [Z/H] ^ +0.3, as found for two bright, inner red GCs as 
well as the integrated galaxy light (Strader et al. 2004; Howell 
et al. 2004). This possibility is illustrated in Figure 9, which 
compares the models shown in Figure 8d with the combi- 
nation of the 1.5 Gyr FZOl model and the Maraston (2005) 
model of 1.5 Gyr and [Z/H] = +0.3. Formally, the least- 
squares fit of this model to the LF indeed yields the lowest 
RMS residuals of all models for intermediate ages considered 
here. For reference. Figure 9 also shows a fit of a 12 Gyr FZOl 
model with a metallicity similar to that of an average Galac- 
tic bulge GC ([Z/H] = -0.7, see Puzia et al. 2002; Maraston 
et al. 2003) to the LF of the inner 50% of the red GCs. As 



TABLE 3 

Fits of FZO 1 mass functions to luminosity function of the 

INNER 50% of the red GCS". 



Age" 


[Z/HY 


rms 


(Gyr) 


(dex) 


error 


1.5 


-1-0.3 


2.38 


1.5 


0.0 


2.43 


3.0 


0.0 


2.48 


6.0 


0.0 


2.57 


12.0 


-0.7 


3.12 



"GCs with V < 26.5 have been considered in the fits. 
"Assumed age of red GCs. 
'^Assumed [Z/H] of red GCs. 



Table 3 shows, this model does not fit the data nearly as well 
as any of the intermediate-age models. Furthermore, if the 
majority of the red GCs would be ~ 12 Gyr old and of mod- 
erate metallicity, one would expect this model to fit the red 
GC LF independent of galactocentric radius, as it does for the 
metal-rich GCs associated with the Galactic bulge as well as 
for red GCs in "normal" giant elhpticals. In contrast, the LF 
of the 50% outermost red GCs is consistent with a power law 
as shown in Figure 8c. 

Another way to test whether the red GC subpopulation is 
more likely to be dominated by an intermediate-age popula- 
tion with [Z/H] > or by an old bulge-like population with 
[Z/H] ^ —0.7 is by comparing the observed mean colors of 
GCs as a function of their luminosity with the predictions of 
the cluster disruption models fit to the LFs of the blue and red 
GCs (i.e., the curves in Figures 8a and 9). This comparison 
is shown in Figure 10 for the inner half of the GC popula- 
tion (r < 46'.'7), for which we showed that the FZOl models 
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Fig. 10. — Mean V — I colors of GC candidates as a function of V mag- 
nitude in the inner 46'.'7. The squai'es represent measured mean V — I colors 
after statistical subtraction of compact background galaxies (see discussion 
in Sect. 4.2.2). The dashed lines indicate mean V — I colors as predicted by 
the FZOl models normalized to the LPs shown in Figures 8 and 9. The line 
types are the same as in those figures, and ages (in Gyr) are indicated on the 
light-hand side of the diagram. Note that the fit of the intermediate-age mod- 
els to the data is significantly better than that of the old model with metallicity 
similar to that of bulge GCs in our Galaxy ([Z/H] = —0.7). 



fit the LFs of both the blue and the red GC subpopulations 
quite well. To calculate mean colors as predicted by the mod- 
els, we assume V — I = 0.93 and V — I = 1.09 for the blue 
and red subpopulations, respectively (cf. Table 2). The ob- 
served mean V ~ I colors are calculated for GC candidates 
with 0.7 < V — I < 1.5. Figure 10 shows that the observed 
behavior of the mean GC color as function of V magnitude 
(the red clusters dominating at both the bright and faint ends 
of the LP) is clearly better fit by metal-rich intermediate-age 
models than by the old model with bulge-GC-like metallicity 
([Z/H] = -0.7). (Incidentally, the 1.5 Gyr model with [Z/H] 
^ +0.3 yields the best fit to the observed mean colors, just as 
it does to the LF of the red GCs.) This result is mainly due to 
the lack of high-luminosity red GCs in the 12 Gyr FZOl model 
with [Z/H] = —0.7 relative to the situation for the models of 
younger age. The underlying reason is that at high GC masses 
and old ages, age fading from one luminosity bin to the next 
(lower) one is somewhat faster than the disruption rate of GCs 
in that lower luminosity (or mass) bin. This is demonstrated in 
Figure 1 1 in which we plot the same model LFs, normalized 
to each other at My ~ —8.5 (in the steep part of the high- 
luminosity wing of the LF), both on logarithmic and linear 
scales. The bright end of the LF moves to lower luminosities 
between 6 and 12 Gyr 

From the above considerations, it seems most likely that 
the red GC subpopulation in NGC 3610 is dominated by an 
intermediate-age population with [Z/H] > 0.0 rather than 
an old population with [Z/H] ^ —0.7. As Figure 9 shows, 
the location of the turnover found for the inner half of the 
red GC system is in rather good agreement with those of the 
FZOl models for intermediate ages and solar metallicity. It 
thus seems reasonable to suggest that the LF of this second- 
generation GC system will evolve further to become consis- 
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Fig. 1 1 . — Predicted V-band GC LFs in relative logarithmic (top panel) 
and linear (bottom panel) units. Then drawn model curves are Fall & Zhang 
(2001) GC disruption models for the following age/metallicity combinations 
(using the Maraston (2005) SSP models): Age = 1.5 Gyr and [Z/H] = 0.0 
(red long-dashed line), Age = 3 Gyr and [Z/H] = 0.0 (magenta short-dashed 
line). Age = 6 Gyr and [Z/H] = 0.0 (black short-dash-dotted line) and Age 
= 12 Gyr and [Z/H] = —0.7 (blue solid line). These age/metallicity combi- 
nations fit the color of the red peak in the GC color distribution fairly well 
and are also plotted in Fig. 9. The curves are normalized to one another at 
My = —8.5, in the steep part of the high-luminosity wing of the LF. Note 
that the luminosity fading predicted by the SSP models of Maraston is some- 
what faster than the GC disruption rate predicted by the Fall & Zhang models 
at high masses. 



tent with the LFs of red GCs in 'normal' giant ellipticals. 

Our results seem to differ from those of Strader et al. (2004) 
who argued (based on spectra taken using the Keck 1 tele- 
scope) that 7 out of their 9 red GCs in NGC 3610 are old 
(> 10 Gyr) rather than of intermediate age. To address this 
apparent discrepancy, we offer the following insights. (;) 
Most red GCs targeted by Strader et al. (2004) were located 
in the outer regions of NGC 3610, whereas the majority of the 
bright red GCs is centrally located, as might be expected for 
GCs formed in a relatively recent starburst event associated 
with a dissipative galaxy merger This is illustrated in Fig- 
ure 12 which compares a histogram of the projected radii of 
the red GCs targeted by Strader et al. (2004) with that of the 
red GC system from the current study, (ii) The faintness of 
the sources and the relatively high galaxy background caused 
many of Strader et al.'s GC spectra to be of marginal qual- 
ity. Applying a commonly assumed absolute minimum re- 
quirement on the signal-to-noise (S/N) ratio to allow a useful 
age/metallicity analysis (S/N > 10 per A, see e.g. Kuntschner 
et al. 2002; Larsen et al. 2003; Beasley et al. 2004; Puzia et al. 
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Fig. 12. — Distribution of projected galactocentric radii of red GCs studied 
by Strader et al. (2004. dashed lines, hatched histogram) compared with that 
of the red GC candidates in the current study (solid lines, open histogram). 
The latter histogram was normalized to the former in the outermost radial 
bin. The dark-hatched histogram (in black) highlights the 4 GCs for which 
Strader et al. (2004) obtained spectra with S/N ratios > 10 per A. 



2002, 2005; Puzia, Kissler-Patig, & Goudfrooij 2006), Strader 
et al.'s results amount to four red GCs with spectra of ade- 
quate S/N of which two are old and two are of intermediate 
age. (This excludes clusters W22 and W28 [i.e., #20 and #26 
in the current study, respectively], which are clearly blue GCs 
as judged from the ACS data, cf. Table 1 .) While we believe 
a small sample of 4 GCs with adequate spectral data does not 
allow one to draw strong conclusions regarding the fraction of 
red GCs being old vs. of intermediate age, it is clear from the 
Strader et al. (2004) study that old GCs do exist among the red 
GC subsample in NGC 3610. As these old, red GCs presum- 
ably stem from the progenitor galaxies that merged to form 
NGC 3610, it seems likely that at least one of those galaxies 
possessed a significant bulge. For example, the Milky Way 
(a giant Sbc galaxy) hosts about 20 - 30 old, metal-rich GCs 
associated with its bulge (e.g., Minniti 1995). The LF ex- 
pected for a population of 50 such GCs according to the FZOl 
model for an age of 12 Gyr is illustrated in Figure 9b as a 
solid curve. Note that the contribution of such a population is 
negligible relative to the observed numbers of red GCs, espe- 
cially since Figure 9c only refers to the inner 50% of the red 
GC subsystem. 

4.3. Radial Distributions and Specific Frequencies 
The specific frequency of GCs, 

(Harris & van den Bergh 1981), i.e., the number of star clus- 
ters per galaxy luminosity in units of My = —15, is an im- 
portant parameter of GC systems in galaxies. Sn is known 
to increase along the Hubble sequence, with mean values 
ranging from {Sn) = 0.5 for Sc spirals (Goudfrooij et al. 
2003; Chandar et al. 2004) to 2.9 for 'normal' giant ellipticals 
outside galaxy clusters, albeit with considerable galaxy-to- 



galaxy scatter (e.g.. Ashman & Zepf 1998; Brodie & Strader 
2006). 

The intermediate-age merger remnant NGC 3610 is a well- 
suited probe to test whether the difference in specific fre- 
quency between spirals and (giant) ellipticals can be ac- 
counted for by secondary populations of clusters created dur- 
ing (gas-rich) mergers. We know that its body has already set- 
tled down to a typical r^/** surface brightness law, we have a 
relatively good estimate of the age of the major merger during 
which the secondary population of GCs were created (both 
from the properties of the GCs themselves and from deep 
spectra of the diffuse galaxy light; Howell et al. 2004), and 
it occurred long enough ago for dynamical evolution of the 
system to have already had a significant impact. Hence, esti- 
mates of evolutionary effects ought to be less uncertain than 
for younger mergers. In the remainder of this section, we de- 
rive the current value of Sn for NGC 3610 and estimate its 
future evolution. 

4.3. 1. Radial Number Density Distribution of Blue and Red 

Clusters 

The radial distribution of the surface number density of blue 
and red GCs in NGC 3610 was derived from the ACS data as 
follows. To avoid incompleteness-related issues, we only con- 
sidered GCs brighter than the 50% completeness limit in the 
central regions {V < 24.5) for this exercise. We calculated 
the completeness-corrected number of blue and red GCs in 
four annuli, logarithmically spaced in projected galactocentric 
radius. Next, we divided the number of clusters in each annu- 
lus by the appropriate area coverage to yield a surface density 
(taking into account the limited azimuthal sky coverage of the 
data at large radii). Finally, we subtracted the surface number 
density of foreground stars and compact background galaxies 
for both blue and red GCs, measured as detailed in Sect. 4. 1 
above. 

Figure 1 3 shows the derived number density distributions of 
the blue and red GCs, along with the surface brightness pro- 
file of the galaxy light for comparison purposes. The surface 
number density profile of the red GCs is largely consistent 
with the galaxy light profile (the innermost surface number 
density value is slightly below the latter), whereas that of the 
blue GCs is significantly flatter than the galaxy light profile. 
Quantitatively, least-square de Vaucouleurs' (1948) r^^^ pro- 
file fits to the GC number density profiles yield — 128 ± 24 
arcsec (or 20.5±3.8kpc) for the blue GCs and re = 19.6±0.2 
arcsec (or 3.14±0.03kpc) for the red GCs. While radial num- 
ber density profiles of red GCs in giant ellipticals are typically 
found to be steeper than those of blue GCs (e.g.. Ashman & 
Zepf 1998), the difference in slope is typically smaller than 
that found here for NGC 3610. In particular, the number den- 
sity profile of the red GCs is significantly closer to that of the 
galaxy light in NGC 3610 than it is in 'normal' giant elliptical 
galaxies (Forbes et al. 1998, 2006; Dirsch et al. 2003, 2005; 
Bassino et al. 2006), where the red GCs are more depleted 
towards the center than the galaxy light. This difference is 
likely due to the cumulative effects of two-body relaxation 
and tidal shocking of GCs (both of which are most effective 
in the central regions as discussed in Sect. 4.2.2 above) and 
dynamical friction (the latter only for the most massive GCs, 
Vesperini 2001). The more moderate flattening of the number 
density of the red GC system of NGC 3610 towards the center 
relative to that of 'normal' giant ellipticals is consistent with 
the intermediate-age nature of NGC 3610. 
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Fig. 13. — Radial number density profiles for the GC system of NGC 3610 
down to y = 24.5. Data have been corrected for completeness, areal cov- 
erage, and contamination by compact background galaxies as explained in 
the text. Filled circles represent the blue GCs while open circles represent 
the red GCs. Dotted lines drawn through the two sets of circles represent de 
Vaucouleurs' r^^^ profile fits to the respective GC number density data. The 
dashed line represents the surface brightness profile of the integrated y-band 
light of the galaxy (taken from Goudfrooij et al. 1994a) shifted by an arbitrary 
constant, while the solid line does the same for the integrated galaxy light as 
measured from the ACS V-band images. 



4.3.2. Specific Frequencies of Blue and Red Clusters 

The number of GCs in the formula for the specific fre- 
quency Sn was originally defined for 'old' GC systems which 
show LFs consistent with a Gaussian, namely as twice the 
number of GCs brighter than the GCLF turn-over magnitude 
(Harris & van den Bergh 1981). This is the procedure we 
follow for the blue GCs, for which the turnover magnitude 
My = —7.2 corresponds to a mass of 2.4 x 10^ Mq accord- 
ing to the Maraston (2005) SSP model for an age of 14 Gyr 
and [Z/H] = -1.7. As the red GCs in NGC 3610 are most 
likely predominantly of intermediate age rather than old, we 
use an alternative definition of Sn for the red GCs: We calcu- 
late twice the number of red GCs more massive than the mass 
associated with the turnover magnitude of the Galactic GC 
system, taking into account the uncertainties of the assumed 
age and metallicity for the red GCs. To this end, we use the 
SSP tables of Maraston (2005) and age/metallicity combina- 
tions ranging from Age =1.5 Gyr and [Z/H] = +0.3 to Age 
= 4 Gyr and [Z/H] = 0.0 (cf. Figure 6). 

Using the r^/^ profile fits to the GC number density profiles 
mentioned above, we calculate the total number of blue and 
red GCs more massive than 2.4 x 10^ Mq. For the purpose 
of Sn calculations, we separate blue and red GCs at V — I = 
1.00 and we assume a distance uncertainty of 15%. Results 
are listed in Table 4. 

Using a total V^-band magnitude of 10.84 for NGC 3610 
(de Vaucouleurs et al. 1991) and hence = -21.81 at the 
assumed distance of NGC 3610, the total present-day Sn = 
1.36 ± 0.58. To allow a proper comparison of this estimate 
with those of 'normal', old ellipticals, this value still needs to 
be corrected for two evolutionary effects: 

1. Luminosity fading of the diffuse light of NGC 3610. 
Treating the diffuse light of NGC 3610 as a SSP, we 
assume an age of 3 ± 1 Gyr for NGC 3610 and an age 
of 10 Gyr for a 'normal' elliptical. This luminosity fad- 
ing amounts to 1.17 ± 0.31 mag in the V band. An 
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Fig. 14. — Histogram of specific frequencies (Sn values) of GCs in ellipti- 
cal galaxies in the compilation of Ashman & Zepf (1998). Central dominant 
galaxies in galaxy clusters have been excluded. The sohd lines represent all 
ellipticals in the compilation, whereas the short-dashed lines represent all el- 
lipticals outside galaxy clusters. For comparison, the specific frequency of 
GCs in NGC 3610 (Sm = 3.8 ± 1.7) is indicated by a vertical long-dashed 
line; a shaded region around it indicates its uncertainty. 



alternative option is that the diffuse light of NGC 3610 
is produced by a combination of an old and a young 
population. Using the Maraston (2005) models at so- 
lar metallicity, we find that the color of the diffuse light 
(B- / = 1.92 and V"- / = 1.11, see Goudfrooij et al. 
1994a) can be modeled as a combination of a 10-Gyr- 
old component (73% by luminosity) plus a 1-Gyr-old 
component (27% by luminosity). In this case, the lumi- 
nosity fading of the young component to an age of 10 
Gyr would amount to 0.80 mag, which is near the low 
end of the range of fading values calculated for the SSP 
age of 3 ± 1 Gyr mentioned above. 

2. The fraction of red GCs more massive than 2.4 x 10^ 
Mq that will be disrupted during this period. To esti- 
mate this fraction, we use the FZOl model mass func- 
tions that involve a power-law initial mass function 
and find that the number of GCs more massive than 
2.4 X 10^ M0 decreases by 12.0% between ages of 1.5 
and 12 Gyr, and by 9.8% between ages of 3 and 12 Gyr 
For the current estimate, we adopt a value of 10%. 

Taking the corrections mentioned above into account, our 
final estimate of the total specific frequency of GCs in 
NGC 3610 at an age of 10 Gyr is Sn ^ 3.8 ± 1.7. The frac- 
tion of red GCs contributing to this value is 0.45±0.05. While 
these values should be checked using a future deep wide-field 
imaging study of NGC 3610 (e.g., the study of NGC 7814 by 
Rhode & Zepf [2003] showed that one can overestimate Sn 
by ^ 20-30% when extrapolating results from HST images 
to large radii), we do note that this value of 5*^ is consistent 
with typical values found for giant ellipticals in loose groups'*. 
This is illustrated by Figure 14 which compares the value for 

* NGC 3610 is the brightest galaxy of a loose group of galaxies (LGG 234. 
5 members; Garcia 1993). 
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TABLE 4 

Results of Specific Frequency calculations. 



Age^ [Z/H]* Nblue'^ Nred"* •S'jv^bluG.now" 'S']V,rcd,now' "Sjv^tot.now* Nred.old* 5iV,bluc,old' Sjv,rod,old^ 'S'Ar,tot,old'' 

1.5 +0.3 376 ±136 412 ± 40 0.71 it 0.36 0.78 ±0.28 1.49 ±0.59 371 ± 36 2.09 ± 1.21 2.06 ± 0.77 4.15 ±1.76 

2.0 0.0 376 ±136 396 ± 40 0.71 ± 0.36 0.75 ± 0.27 1.46 ±0.58 356 ± 36 2.09 ±1.21 2.02 ± 0.75 4.11 ± 1.75 

3.0 0.0 376 ±136 304 ± 36 0.71 ± 0.36 0.57 ±0.21 1.28 ±0.52 274 ±32 2.09 ± 1.21 1.52 ± 0.58 3.61 ± 1.61 

4.0 0.0 376 ±136 268 ± 32 0.71 ± 0.36 0.51 ±0.19 1.22 ±0.51 241 ± 29 2.09 ± 1.21 1.34 ±0.51 3.43 ± 1.56 



"Assumed age of red GCs in Gyr. 

''Assumed [Z/H] of red GCs in dex. 

'^Total number of blue GCs around NGC 3610. 

''Total number of red GCs around NGC 3610. 

■^Present-day specific frequency of blue GCs. 

^Present-day specific frequency of red GCs. 

STotal present-day specific frequency of GCs. 

''Number of red GCs at age of 10 Gyr, corrected for disruption. 

'Specific frequency of blue GCs at age of 10 Gyr. 

J Specific frequency of red GCs at age of 10 Gyr. 

•'Total specific frequency of GCs at age of 10 Gyr. 



NGC 3610 found here with a histogram of Sn values for 
'normal' elUpticals (taken from the compilation of Ashman 
& Zepf 1998, which does not include NGC 3610). Quantita- 
tively, the mean and standard deviation of specific frequencies 
of ellipticals outside galaxy clusters in the latter compilation 
is Sjv = 2.9 ± 1.9. 

5. SUMMARY AND CONCLUSIONS 

We have presented a study of globular cluster candidates 
in NGC 3610, a giant elliptical galaxy that is very likely 
an intermediate-age merger remnant, using deep BVI images 
taken with the ACS camera on board the Hubble Space Tele- 
scope. Our main results are as follows. 

1 . The globular clusters show a clear bimodal color dis- 
tribution, both in i? — / and in V — I. The peak colors 
of the blue subpopulation are B — / = 1.61 ± O.Oland 
V — I = 0.93 ± 0.01, and those of the red subpopula- 
tion are B-/ = 1.92 ± 0.01 and V-I = 1.09 ± 0.01. 
The colors of the blue clusters are consistent with those 
of blue clusters in other elliptical galaxies as well as 
of halo clusters in our Galaxy, while the colors of the 
red clusters are significantly bluer than the colors in- 
dicated by the relation of red cluster color with galaxy 
luminosity among 'normal' elliptical galaxies. Given 
the additional fact that the brightest red clusters are up 
to 0.5 mag brighter than the brightest blue clusters in 
NGC 3610, we interpret this as evidence for a younger 
age of the red clusters in NGC 3610 relative to normal 
giant elhpticals rather than a significantly lower metal- 
licity. 

2. The luminosity function of the blue, metal-poor sub- 
population of clusters in NGC 3610 is similar to a Gaus- 
sian with ~ —7.2 and a ~ 1.2, consistent with 
values typical for old blue cluster populations in nor- 
mal elliptical galaxies. 

3. The luminosity function of the full red subpopulation 
of clusters in NGC 3610 is consistent with a power law 
(l){L)dL oc L"dL with a = -1.78 ± 0.06, which is 
consistent with the result of W02 using earlier WFPC2 
data. However, the luminosity function of the inner 



half of the red cluster subsystem shows a clear flat- 
tening consistent with a turnover at My ~ —6.2. A 
comparison with predictions of GC disruption models 
in conjunction with simple stellar population models at 
various ages and metallicities shows that both the lumi- 
nosity function of the inner half of the red cluster sub- 
system and the relation of mean cluster color with lumi- 
nosity are best fit by a population of metal-rich clusters 
of intermediate age (~ 1.5 — 3 Gyr). As similar results 
were found for the red GCs in the 3-Gyr-old merger 
remnant NGC 1316 (Goudfrooij et al. 2004), this re- 
sult supports the assertion that dynamical evolution of 
metal-rich GC populations formed in gas-rich galaxy 
mergers indeed changes their LF properties to become 
consistent with those of the red, metal-rich GC subsys- 
tems that are ubiquitous in 'normal' giant ellipticals. 

4. The surface number density profile of the blue GC sub- 
population is much more radially extended than the sur- 
face brightness profile of the galaxy hght, similar to the 
situation in 'normal', old elliptical galaxies. In con- 
trast, the radial profile of the red GC subpopulation fol- 
lows that of the galaxy light almost all the way into the 
galaxy center. The difference in steepness between the 
radial profile of the red GC population and that of the 
galaxy light is smaller than reported for 'normal', old 
elliptical galaxies. The combination of this result with 
the finding that the red GCs and the galaxy light ap- 
pear to have ages and metallicities consistent with each 
other suggests that most of the red GCs were formed 
during the same event as the bulk of the field stars in 
NGC 3610 and that they have suffered Uttle disruption 
after violent relaxation with respect to the situation in 
'normal' ellipticals. 

5. We determine the specific frequency of clusters in 
NGC 3610 and find a present-day value of Sat = 1.4 ± 
0.6. Using pubhshed age estimates for the diffuse light 
of NGC 3610 as well as cluster disruption models, we 
estimate that this value will increase to Sn = 3.8 ± 1.7 
at an age of 10 Gyr, which is consistent with typical ^jv 
values of 'normal' elUptical galaxies in loose groups. 



HST/ACS Imaging of Globular Clusters in NGC 3610 
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Our results support a picture in which the formation process 
of giant ellipticals with significant populations of metal-rich 
GCs was similar to that of gas-rich galaxy mergers observed 
today. 
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